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Abstract 

A systematic method to account for electron correlation in periodic systems which 
can predict quantitatively correct band structures of non-conducting solids from first 
principles is presented. Using localized Hartree-Fock orbitals (both occupied and 
virtual ones), an effective Hamiltonian is built up whose matrix elements can easily 
be transferred from finite to infinite systems. To describe the correlation effects 
wave-function-based multireference configuration interaction (MRCI) calculations 
with singly and doubly excited configurations are performed. This way it is possible 
to generate, both, valence and conduction bands with all correlation effects taken 
into account, irans-polyacetylene is chosen as a test system. 
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1 Introduction 

Recently, with the rapid progress in computing facilities, rather sophisticated 
numerical methods started to arise which attempt to determine quantitatively 
correct electronic band structures of solids from first principles. These methods 
can be grouped in two families: those focusing on ground-state properties such 
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as the electron density (or the quasi-particle Green function) and those trying 
to approach the full many-body electron wave function. 

As a rule, traditional density functional methods are used in solid state physics 
to study the electronic properties of infinite periodic systems, but since these 
approaches are not designed for excited states they usually strongly underes- 
timate the band gaps of non-conducting systems. Also these methods do in 
general not allow improvements in some controlled way. 

On the other hand, methods based on many-body wave functions of finite 
systems are well-established in quantum chemistry. They usually start from 
an approximate wave function (e.g. the Hartree-Fock (HF) wave function) and 
approach the correct solution by improving the wave function systematically. 
The development of such methods for solids only started recently (see below). 
The band gaps obtained on the HF level (via Koopmans' theorem) are usually 
far too large, but the situation improves steadily upon successive inclusion of 
correlation effects. In this sense the wave-function-based methods are more 
general though more demanding in predicting quantitatively correct electronic 
structures of periodic system than the common density functional methods. 

Among the recently developed wave-function-based approaches to correlated 
band structures one can distinguish those which employ local Hamiltonian ma- 
trix elements in real space [1,2,3,4], those which work directly in A:-space and 
use many-body perturbation theory to account for electron correlation [5,6], 
and those which use an approximate Green function formalism [4,7,8]. Other 
approaches to electron correlation in solids [9,10] are presently implemented, 
corroborating the currently high interest in this new and rapidly growing re- 
search field. 

The aim of this Letter is to present the local Hamiltonian method which 
originates from the approach to valence bands described in Ref. [1,2,3] but 
which can now also treat conduction bands and thus allows to predict com- 
plete band structures (and in particular band gaps) of non-conducting systems 
with controlled accuracy. The presented method is designed to be used with 
any suitable standard quantum-chemical program package. We decided for 
the multireference version of the single and double configuration interaction 
method (MRCI(SD)) to account for electron correlation in both the ground 
and excited states of the system together with a special size-consistency cor- 
rection we have developed for that purpose [11]. This explicitly includes up 
to five-particle configurations (three holes and two electrons or vice versa) 
and therefore goes far beyond any of the other wave-function-based methods 
mentioned above. 

To achieve our goal, the original problem of determining delocalized electronic 
states with a crystal momentum (Bloch states) is reformulated in terms of 
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localized orbitals in real space (Wannier functions). In this representation one 
can take advantage of the predominantly local character of correlation effects 
and evaluate them step-by-step in finite clusters of the periodic system. The 
approach is purely ab initio. We adopt the experimentally known geometry of 
our test system (£rans-polyacetylene) and we do not make any approximation 
for the electron-electron or electron-nuclear interaction. 



2 The method 



The valence band energies e kv of a periodic system are defined as 

—e kv = IP ku = E^' 1 - Eq (1) 

where Eq is the ground-state energy of the neutral iV-electron system and 
E*- 1 = {^ua l \H\^kua) is the energy of that excited state V^' 1 of the 
(iV-l)-electron system which corresponds to a quasi-hole of spin —a, crystal 
momentum — k and valence band index v. The corresponding equation for the 
conduction band energies e kfJi reads 

-£ fcM = EA fcM = Eq - E^ 1 (2) 

with \x labeling the conduction bands and E^ 1 = (^fc^l^l^fc^ 1 ) referring 
to the excited-state of the (iV+i)-electron system which is associated with a 
quasi-particle of spin a in a virtual Bloch orbital f kfl . Note, that the band 
energies are differences of total energies rather than simple one-particle en- 
ergies. Also note, that excited states which do not have any counter part in 
the single-orbital picture such as satellite and shake-up states are explicitly 
excluded here. 

On the Hartree-Fock level the ground-state wave function of the neutral system 
is approximated by a single Slater determinant $, and because in infinite 
systems there is no orbital relaxation upon removal or adjunct of a single 
electron, Koopmans' theorem holds such that the excited states of the (iV±l)- 
electron system are simply given by 

K~ 1 )=c k ^m and K + ;) = d$) . (3) 

Here c kva and c k create or annihilate an electron in the respective HF spin 
orbitals which can directly be obtained (together with the HF band energies) 
from standard HF program packages for periodic systems like CRYSTAL [12]. 

Applying separate multi-band Wannier transformations to the occupied and 
virtual Bloch orbitals provides two sets of local Wannier functions, ipun and 
^Rmi where R denotes the lattice vector to the unit cell associated with each 
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Wannier function (see Fig. 1). Using these Wannier functions so-called local 
one-particle configurations 

\*r£) = CRn*\*) and |*K> = 4j*> (4) 

can be introduced, where CR ncr and c Rma create or annihilate an electron from 
the respective Wannier spin orbital. Having switched to the local configura- 
tions the Hartree-Fock band energies can be recovered by diagonalizing the 
/c-dependent 'ionization potential' and 'electron affinity' matrices, 

IPnn F '(fc)=E eifcRlP 3(^) With 
R 

IPZ(R) = (KnMl^a) ~ E» 5 OR 8 nn , (5) 

and 



EA« F m , (fc) = £ e ikR EAZAR) with 

R 

EAZ ll (R)=E£6 OR 6 mm , - (K^\H\$r+1) • (6) 



The same holds for the correlated band energies provided the Wannier trans- 
formations which mediates between the Bloch and the Wannier orbitals is 
applied to the excited-state wave functions as well, yielding 

C^WEe"* 8 *^ and (7) 

C=EWEe M *Si . (8) 

m R 

where U vn {k) and U^ m {k) are the band- mixing matrices of the two Wannier 
transformations. The correlated ionization potential and electron affinity ma- 
trices IP nn ' and EA mm ' are defined in perfect analogy to Eqs. (5) and (6). 

The transformed many-body wave functions and VP^* are sometimes 

referred to as 'correlated local holes' and 'correlated local electrons'. Yet, like 
the Wannier functions or any other kind of localized orbitals they are pure 
mathematical constructions only meant to provide a convenient local repre- 
sentation of the delocalized eigenstates ^^~ a l and ^j^ 1 one ls looking for. 

The benefit of this local representation is two-fold. Firstly, the off-diagonal 
(or 'hopping') matrix elements of lP nn >(R) and EA mm /(i2) decay rapidly with 
increasing distance R when well-localized Wannier orbitals are used such that 
the (in fact infinite) summation over R in Eqs. (5) and (6) can be truncated at 
some distance R cu t- Secondly, the individual local matrix elements lP nn i(R) 
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Fig. 1. Two typical local orbitals in irans-polyacetylene. Contours are in steps of 
0.025 au. 

and EA mm /(i2) can be evaluated in finite clusters of the periodic system, 
because excitations from the involved local spin orbitals into some distant 
spin orbitals become negligible once they are sufficiently apart. In addition, 
as outlines in Ref. [2,3], asymptotic corrections can be applied to account for 
the missing contributions, if necessary. 

To obtain the local matrix elements (LMEs) from the clusters our local Hamil- 
tonian approach is applied to the finite cluster as well leading to a set of local- 
ized occupied and virtual orbitals ip a and (p r and the corresponding molecular 
ionization potential and electron affinity matrices IP aa ' and EA rr /. However, 
instead of determining the LMEs in order to find the band energies as done 
for the periodic system, one now proceeds in the opposite direction. Starting 
from the excited-state energies E+ and E~ of the cationic and anionic cluster, 
which can be calculated (together with the molecular ground-state energy E ) 
by means of any suitable quantum chemical correlation method, the unknown 
LMEs are evaluated according to 

IP*,, = J2 K U„aU: a , - E 5 aa , and (9) 
EA^ = E 6 rr ,-'EE-U^U lir , ■ (10) 

The unitary matrices U va and U^ r stem from some suitable localization such as 
the Foster-Boys procedure [13] (see Fig. 1) and are the molecular counter parts 
of the band-mixing matrices. We used the MRCI(SD) option [14,15,16] of the 
MOLPRO program package [17] taking the unrelaxed Hartree-Fock orbitals of 
the neutral system to expand the wave functions of the charged molecules as 
well. The lack of size-consistency (SC) inherent to truncated CI schemes could 
be overcome by developing a special SC correction for the (iV±l)-electron 
systems [11] similar to the well-established Pople correction [18] for closed- 
shell systems. To facilitate the still quite demanding cluster calculations the 
incremental scheme was employed (see Ref. [2,3,11] for details). 

Despite the formal hole-particle dualism in the equations above, there is quite 
a difference between valence and conduction band calculations for ab initio 
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Hamiltonians. Virtual orbitals are much harder to localize than occupied ones 
and the anionic clusters tend to become electronically unstable such that some 
embedding or confinement techniques have to be applied to make the calcu- 
lations feasible. 



3 Results 

As test system for our method we chose a irans-polyacetylene single chain, an 
infinite, flat hydro-carbon molecule with alternating C-C bond lengths and 
C2H2 as repeat unit (see Fig. 1). The experimental geometry of the carbon 
skeleton [19], rf(C-C) = 1.45A and d(C—C) = 1.36A with a lattice constant 
of 2.455A, is adopted resulting in a bond angle Z(C-C,C=C) of 121.74°, and 
the C-H bond length is 1.087A (taken from Ref. [20]). Flexible correlation- 
consistent valence triple zeta (VTZ) basis sets [21] were used throughout, to 
allow reasonable correlation calculations also for the anionic (A^+l)-electron 
states. 

X-ray scattering tells [22] that crystalline trans-polyacetylene (tPA) consists 
of weakly coupled fishbone-like packed tPA single chains all aligned in the 
same direction, as is also evident from the high anisotropy of the electronic 
properties of bulk tPA. Our HF calculations with CRYSTAL give valence band 
widths of 5-10 eV for crystal momenta k parallel to the chains but not more 
than 0.4 eV in the perpendicular directions. The same holds for the low-lying 
conduction bands up to about 3 eV above the lowest conduction band energy. 
Hence, up to this energy, bulk tPA can, in fact, be considered as a quasi one- 
dimensional crystal and a single tPA chain with the finite basis sets is an 
appropriate model for investigating that part of the electronic structure. 

The band structure of a single tPA chain as obtained by our local Hamilto- 
nian approach is depicted in Fig. 2 where it is compared to the HF data. All 
five valence bands and the three lowest conduction bands of tPA are shown. 
Evidently, the band structure changes dramatically upon inclusion of electron 
correlation. The valence bands shift upwards, the conduction bands down- 
wards, and a flattening of all bands (significantly more pronounced in the 
valence region) is observed leading to the desired narrowing of the band gap. 
Other interesting but more subtle changes are also discernible, e.g. the lifting 
of accidental degeneracies or the narrowing of avoided crossings. The band 
gap of tPA reduces by 2.31 eV from 6.42 eV on the HF level down to 4.11 eV 
when correlation is included. This value is perfect in line with similar values, 
2.35 and 2.38 eV, reported in two independent studies on tPA single chains 
[5] and [6], where the tt bands alone were investigated by means of M0ller- 
Plesset perturbation theory (MP2). The lowest ionization potential becomes 
4.58 eV, the highest electron affinity 0.47 eV, a positive value finally, indicat- 
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Fig. 2. Correlated (solid line) and Hartree-Fock (dashed line) band structure of an 
infinite irans-polyacetylene single chain. 

ing that trans-polyacetylene is able to retain an extra electron as is known 
from experiment [23]. 

The accuracy of our approach is governed by the error in the determination 
of the LMEs by the incremental scheme. For tPA it mainly arises from the 
truncation of the incremental series after the second order corrections. The 
discarded long range polarization effects turned out to be negligible. The re- 
sulting uncertainty in the band structure could be estimated to not exceed 
0.3 eV for any Appoint which is less than 1% of the total band energy range. 
A more detailed error analysis can be found elsewhere [11]. The new size- 
consistency (SC) correction for the excited states is crucial here to reduce 
the SC error in the LME increments to 0.01% of the contributing correlation 
energies. The truncation error of the lattice summations in Eqs. (5) and (6) 
was analyzed on the HF level. Keeping only LMEs above 1 mHartree (which 
results in 59 LMEs for the valence and 39 for the conduction bands) the band 
energies from a canonical HF calculation with CRYSTAL could be reproduced 
within 0.2 eV. These deviations can easily be accounted for by adding them 
to the band energies emerging from the LMEs (as done in Fig. 2). 



4 Discussion and Conclusions 



In optical studies on bulk trans-polyacetylene [22,24,25] an absorption max- 
imum at about 2 eV is observed. Taking into account the evidence, both 
experimentally and theoretically for the existence of an excitonic gap state 
[5,26,27,28], 0.4 eV below the conduction band [27], the fundamental band 
gap of tPA becomes 2.4 eV, which is noticeably smaller than the value of 4.1 
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eV we (and others [5,6]) found. 7r-band splitting due to interchain coupling was 
made responsible for that discrepancy [6], but for the experimentally found 
fishbone-like stacking of bulk tPA we only find a reduction of the fundamental 
band gap of 0.6 eV on the HF level and suppose (based on the MP2 data in 
Ref. [6]) a further 0.2 eV reduction upon inclusion of correlation leading to a 
gap of 3.3 eV. The more pronounced reductions (~ 1.4 eV) reported in Ref. [6] 
for two tPA chains are considered to be an artifact of the vertical tPA packing 
assumed there. 

The missing long-range interchain polarization existing in the bulk material 
could be the reason why wave-function-based methods applied to tPA single 
chains (or stacks of a few of them) consistently overestimate the fundamen- 
tal band gap of tPA, but it could also be that effects beyond bare electron 
correlation have to be included to achieve accordance with experiment. 

Nevertheless, this study demonstrates that our local Hamiltonian approach 
provides a sought scheme to perform reliable wave-function-based correlation 
calculations for excited electron hole and attached-electron states in solids and 
polymers with any suitable quantum chemical correlation method available, 
as corroborated by similar studies in our group on diamond and hydrogen 
fluoride chains. 
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